Long lived acoustic vibrational modes of an embedded nanoparticle 
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Classical continuum elastic calculations show that the acoustic vibrational modes of an embedded 
nanoparticle can be lightly damped even when the longitudinal plane wave acoustic impedances 
Zo = pVL of the nanoparticle and the matrix are the same. It is not necessary for the matrix 
to be less dense or softer than the nanoparticle in order to have long lived vibrational modes. A 
corrected formula for acoustic impedance is provided for the case of longitudinal spherical waves. 
Continuum boundary conditions do not always accurately reflect the microscropic nature of the 
interface between nanoparticle and matrix, and a multi-layer model of the interface reveals the 
possibility of additional reduction of mode damping. 

PACS numbers: 62.20.-x,43.20.Ks,62.25.+g 
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A classical vibrational mode weakly coupled to a cold 
environment will have an amplitude exponentially decay- 
ing with time. The vibrational frequency of the mode 
may be slightly shifted from what it would have been 
if the environmental coupling is reduced to zero. When 
coupling is stronger, the mode loses its individual iden- 
tity, and it becomes necessary instead to consider normal 
modes of the entire system. 

Nearly spherical clusters of atoms a few nanometers 
in size (nanoparticles) have mechanical vibrations which 
can be initiated through sudden thermal excitation with 
a brief laser pulse. [l| Subsequent evolution of the vibra- 
tion can be followed by observing the amount of light 
absorption by a probe pulse following the pump pulse by 
a controlled time. 

It is rarely possible to perform such experiments on iso- 
lated nanoparticles. Normally there is mechanical cou- 
pling to an environment. This can be (1) a powder of 
nanoparticles (where adsorption is important) (2) a sub- 
strate to which the nanoparticle is attached (3) liquid 
solution or (4) a solid surrounding matrix, usually glass. 
Propagation of acoustic waves away from the nanoparti- 
cle is often the dominant mechanism for mechanical en- 
ergy loss. As a result, vibrational modes are damped. 

The organization of this paper is as follows: We begin 
with the current "state of the art" of vibration eigen- 
mode calculation for free and matrix-embedded spheres. 
Then we focus on a particular material combination 
which gives rise to very low damping of the fundamen- 
tal "breathing" mode even though the longitudinal plane 
wave acoustic impedances in the two materials are the 
same. We then derive the correct expression for the 
acoustic impedance for such modes and explain this re- 
sult. Finally, we extend the calculation to more realistic 
materials where the sphere-matrix contact is not perfect 
and to modes having different symetries. 

A quantum mechanical analysis of energy loss of a 



nanoparticle weakly coupled to a substrate heat bath has 
been done.Q Otherwise, the nanoparticle has been ap- 
proximated as a classical isotropic homogeneous linear 
elastic continuum object in contact with a similarly ideal- 
ized zero temperature matrix. We consider such a spher- 
ical nanoparticle of radius Rp, density pp and speeds of 
sound VLp and VTp embedded within a macroscopic ma- 
trix of density p,„ and speeds of sound VLm and VTm ■ Ide- 
alized continuum boundary conditions |3| of continuous 
displacement and normal stress components have been 
applied where the nanoparticle contacts the matrix, but 
these will be reaxamined later on. 

"Pseudomodes" of a nanoparticle embedded in a solid 
or fluid matrix can be found 4] yielding complex val- 
ued frequencies w^^™ of the "complex frequency model" 
(CFM). q can be torsional (TOR) or spheroidal (SPH). 
£ is the angular momentum and m is the z-component. 
n > is the mode index. Re{uj^[^) = 2tt/T are shifted 
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of a free sphere model (FSM) as first calcu- 



lated in 1882 by Lamb^S] and some new modes appear. |y| 
Im{uj^f^) = 1/t = ■^Aoj where r is the damping time 
and Au! is the full width at half maximum (FWHM) of 
a Brillouin or Raman peak. The quality factor of each 
pseudomode is defined as Q = lu/Alu. This paper is 
mainly focused on (SPH,£ = 0) modes which are the 
only ones to have been experimentally observed in pump- 
probe experiments. These modes are also Raman active. 

Interpretation of CFM is not straightforward, since 
the resulting pseudomodes are not orthonormalizable and 
blow up exponentially with the radial coordinate. This 
makes a correspondence to quantum theory unclear. An 
alternative is to calculate all of the normal modes of 
a nanoparticle ("core") surrounded by a large spher- 
ical matrix ( "shell" ) , imposing an arbitrary condition 
on the matrix outer surface, [if The resulting core-shell 
model (CSM) is conservative and has real valued mode 
frequencies. |y| However, there is a continuum of normal 
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modes in the macroscopic limit that the matrix radius is 
infinite. Note that highly damped "matrix modes" are 
not considered in this work. 

The connection between the continuum of CSM modes 
and the discrete spectra of CFM pseudomodes is clearly 
seen by making a plot of the mean squared displace- 
ment within the nanoparticle {u'^)p as a function of fre- 
quency for CSM modes. The resulting peaks corre- 
spond closely to i?e(w^^™) and the FWHM of the CSM 
peaks corresponds closely to 2Im{uj^™). 

Even without applying the complexities of CFM or 
CSM, the qualitative effect of the matrix in shifting and 
broadening FSM frequencies can often be correctly antic- 
ipated by comparing the longitudinal plane wave acoustic 
impedances (LPWAI) Zop = PpVLp and Zom = PmVLm- 

The "rule of thumb" is that nanoparticle vibrational 
modes will be lightly damped when the LPWAI's satisfy 
Zop ^ Zom, or Zop ^ Zom- Otherwise, similar LPWAI's 
of the nanoparticle and matrix allow vibrational energy 
to freely flow out, resulting in strong damping. 

Here we report dramatic exceptions to this. It is pos- 
sible for a nanoparticle to have lightly damped normal 
modes without LPWAI's Zom and Zop being very differ- 
ent, with the extreme possibility of r being large even 
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Furthermore, lightly damped nanoparti- 



cle frequencies can differ greatly from w^^^. These pos- 
sibilities have remained obscure up to now because ex- 
periments were limited primarily to situations in which 
the LPWAI's satisfied Zop > Zom for which FSM works 
extremely well. Furthermore, nanoparticle size variation 
within a sample creates extrinsic broadening of the modes 
which can sometimes exceed intrinsic broadening. 

An example will illustrate when such unexpected re- 
sults can occur: a gold nanoparticle embedded in a di- 
amond matrix. We use the elastic constants of bulk 
fee gold at 300 K0 and use directional averaging^ 
to obtain longitudinal and transverse speeds of sound 
of 3330 m/s and 1250 m/s respectively. We use elas- 
tic constants for bulk crystalline diamond at 300 Kpllj 
and similarly obtain speeds of sound 18190 m/s and 
12350 m/s. The resulting LPWAI's are identical: Zop « 
Zom ~ 6.4 • lO^kg m^^s^^. In the absence of any ap- 
parent impedance mismatch, the "rule of thumb" says 
that the vibrational modes of the nanoparticle would be- 
come so broadened by the presence of the matrix as to 
be unrecognizable. 

Figure n shows something completely different. (u^}p 
for CSM (SPH,£ — 0) modes is plotted versus frequency. 
For comparison, lo^^ are plotted showing Re{uj) as a 
vertical arrow with horizontal arrows to indicate Im{ijj) 
which is half the FWHM. The resulting CSM peaks and 
i?e(ijj^™) match very closely. The gold nanoparticle has 
several lightly damped modes of vibration. 

Figure El shows how misleading naive comparison of 
LPWAI's can be. It was obtained by varying both pm and 
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FIG. 1: Mean square displacement within the nanoparticle 
interior of CSM (SPH,^ = 0) vibrational modes for Au (radius 
i?p=10 nm) in diamond (radius Rm ~ 2000 
indicate positions and FWHM obtained with the CFM. 



only when v^m < I'Lp- Note that for a gold nanoparticle 
in diamond, VTm/vLm ~ 0.68 and Pm/pp ~ 0.18. For 
low matrix mass density, the frequency approaches that 
of the bound sphere model (BSM), ti.'gpHoo' which is for 
a sphere with zero displacement boundary conditions. 

For a purely longitudinal spherical wave in an isotropic 
elastic medium, the acoustic impedance is Z — —Urr/vr 
where stress component Urr — 2pv^err + pi^L — 



this case strain e,- 
radial displacement is Ur 



t,) and velocity Vr — dur/dt. For 
dur/dr and eee — e,f,^ — Ur/r. The 
-'"*(d/dr)(e^'=^7r) where 



vi^k^. The resulting expression for impedance is 



Z{s) 



p{v\ + 4w|(is - s^)) 



wl(1 -I- is) 



(1) 



where s = l/(fcLr). This can only be applied to the 
case of (SPH,^ = 0) modes. In the plane wave limit 
(large r) Z ^ pvL- For a fluid medium, vt = so that 
Z = pvl/O- + is), a result familiar for sound waves (llj. 
VLm while keeping Zom = Zop- Mode broadening is large Another interesting case is Poisson ratio v = ^ (approx- 
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FIG. 2; Normalized CFM frequency (line, left axis) and 
1/Q = Acj/tj (dots, right axis) of the fundamental (SPH,^ = 
0) mode along Zom/Zop = pmVLm/ ppVLp = 1 where 
VTp/vLp = 0.375 (gold) and VTm/vLm = 0.6. Acj is the full 
width at half maximum. 

imately true of many materials) where vt — so that 
Z = pVL to linear order in s. Note that vt /vl = \J \l^2v • 

If s is small, then Re{Z) « pv^ and Im{Z) « 
spvl{'^{vt/vlY — !)• Note that if vt ~ \vl then Im{Z) 
is approximately zero. 

Define a — VLm/vLp- Let Zp(l/^) and Z,„(q;/^) denote 
the acoustic impedance evaluated at Rp of the nanoparti- 
cle and matrix respectively for outgoing waves of dimen- 
sionless frequency ^ = ujRp/vLp- Note that for a given 
pseudomode, these impedances do not depend on Rp. If 
Zp(l/^) = Zjn{a/^) then a longitudinal travelling wave 
can cross the nanoparticle-matrix interface without any 
reflection. 

Consider Zp{l/^): for the (SPH,£ = 0,n) FSM mode, 
s(r — Rp) ^ (n+i)TT - ^'^^ fundamental (n = 0) mode 
s{r — Rp) is not small and Zp{l/£,) is complex- valued. 
But for higher order modes s becomes small so that 
~^ PpVLp- Thus, it is expected that the acoustic 
impedance mismatch at the nanoparticle-matrix interface 
will be different for higher order modes. Figure [T] con- 
firms this. Also, as predicted in a previous paper the 
area under each pseudomode is roughly constant in the 
{u^)p plot. This is important because it means that less 
damped modes compensate their lower frequency broad- 
ening with a higher mean square displacement and there- 
fore can have similar electron-phonon coupling inside the 
nanoparticle as more damped modes. 

Consider the longitudinal radial (r) velocity within the 
nanoparticle {p) to be the sum of inward (i) and outward 
(o) travelling waves: 

vrpoix) = e-*"*Ap„^e"((l/x) + {i/x^)) (2) 
vrp^{y) = e~'^'Ap,u;e'y{{l/y) + {i/y')) (3) 



where x = uir/v^p and y = —ujr/vLp. Apo and Api are 
amplitudes. To avoid singularity at the origin, Apo = 
—Api. Matching of boundary conditions at the interface 
at Rp corresponds to the equation: 

Zp{l/^) - Zmja/O ^ Vrp^i-0 

Zp 

Vrpo 

The roots of this equation are the complex valued dimen- 
sionless frequencies ^^p]^„. This equation is equivalent 
to previous solutions 4] of CFM for the (SPII,0) case, but 
clearly illustrates the role of the acoustic impedances in 
determining the mode frequencies. 

Previous studies[l^ were restricted to analyzing the 
situation of a soft, light matrix where < v^p and 
Pm < Pp, and found that t does not depend on n. But 
we find that in other situations damping can vary dra- 
matically with mode order, as clearly seen in Fig. ^ 

In the matrix, s = a/^ can be small if VLm VLp- 
But if ULm ^ VLp then s can be large, and in this limit 
Zm{ct/C) ^ ^iPm.v'^tml i^Rp), and the effect of the matrix 
is a purely restorative force, as if the matrix is a massless 
spring. Energy cannot be absorbed by the matrix in this 
limit, so that r — > oo when VLp/vLm 0. 

We have explored the lowest (SPII,£ = 0) mode (ig- 
noring highly damped "matrix modes" when they are 
present) in the parameter space of Pp, v^p, vxp, ^Lm 
and and found some general patterns. Notably, good 
impedance matching between nanoparticle and matrix as 
exhibited in very strong lowest (SPH,0) mode damping is 
an unusual occurence requiring special conditions. Close- 
ness of the LPWAI's alone is not sufficient. 

If Vm <^ Vp (where u,n and Vp arc the Poisson ra- 
tios of the matrix and particle respectively) then there 
is never strong damping, but we can have strong damp- 
ing if Vm S> Vp such as for a liquid matrix. Assuming 
the Poisson ratios are close, it is usually necessary that 
Pm > Pp and VLm < VLp- Under these conditions, s is 
small in the matrix and the LPWAI's are close to the true 
impedances. In addition, if Vm and Vp are both close to 
■| then strong damping is still possible if pm is as little 
as Q.Spp. Conversely, light damping is always seen when 
any of the following conditions is satisfied: (1) Vp ^ Vm 

(2) p,n < 0.8pp (3) VLm > l-^VLp. 

As was shown in a previous work, J] the microscopic 
details of the interface {e.g. the specific nature of the 
mechanical connection of the gold crystal lattice to the 
diamond matrix) can tend to isolate (i.e. decouple) the 
nanoparticle from the surrounding medium. In particu- 
lar, continuum-type boundary conditions applied to the 
interface may not be optimal. One rough-and-ready 
way to adapt the CFM and CSM models to tentatively 
take this into account is to introduce an intermediate 
shell layer (X-layer) of inner radius Rp — (ix/2, outer ra- 
dius Rp + (ix/2, density px and speeds of sound vlx 
and vtx- Depending on the elaboration process, the 
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nanoparticle and the surrounding matrix do not neces- 
sarily tightly adhere. Different expansion coefficients for 
example could result in a void in the matrix bigger than 
the nanoparticle and the nanoparticle being only par- 
tially in contact with the matrix. The X-layer is the 
most elementary way to explore the effect of deviations 
from continuum boundary conditions. Such calculations 
are relatively easy to perform using CFM or CSM. 

It is important to demonstrate the robustness of CFM 
results against the microscopic details of the interface. 
The nature of the interface is unknown, but it can be 
guessed to be of lower density, poorly-ordered and not 
strongly-bonded. It is plausible to assume a low density 
and low elastic constants comparable to water or plastic. 
If dx is thin, neither px nor vtx have a strong effect. 
What is important is the elastic constant of the X-layer 
in terms of pxv\x relative to the thickness of the X- 
layer, dx- It is convenient to introduce the dimensionless 
parameter ax = dxPpVLp/{RpPxvj^x)- For examples we 
looked at, the larger ax is, the larger r and Q become. 

A possible example of this may have been seen for gold 
nanoparticles embedded in Ti02. This is unrelated to our 
previous example of Au in diamond. Experimental ob- 
servations of nanoparticle "breathing modes" (those with 
purely radial vibration) for Au nanoparticles in a TiO? 
matrix using the femtosecond pump probe technique [fSj 
reveal surprisingly light damping. Qian et al. obtained 
an oscillation frequency of 0.15 THz (5 cm~^) with a 
damping time of 55 ps {Q ~ 26) for Rp = 10 nm. We use 
Ti02 sound speeds of 8610 m/s and 5160 m/s and density 
4.097 g/cm3.[lj Our XCFM calculations for very small 
ax predict r « (10.7-h21.0crx) ps or Q w 5.92 -h ll.TtJx. 
This T provides only an upper bound on the experimental 
damping time because of dephasing due to particle size 
variation. The large t of Au in Ti02 could be obtained 
by a water-like layer between the nanoparticle and the 
matrix on the order of 0.1 nm thick. Thus, details of the 
nanoparticle-matrix interface at the atomic level are cru- 
cial in determining the damping of embedded nanopar- 
ticle vibrational modes. However, proper first-principles 
quantitative modelling of the interface requires consider- 
ation of detailed features of crystal and electronic prop- 
erties at the atomic level, in particular the nature of the 
bonding. ^3 For the systems we consider, there is insuf- 
ficient information about the interface to attempt such a 
calculation, especially since the matrix is amorphous. 

The ability to choose the matrix so as to have nanopar- 
ticles with lightly damped vibrational modes is impor- 
tant for several reasons. First, by reducing homogeneous 
broadening, it allows inhomogeneous broadening effects 
to stand out clearly. This results in a clearer picture for 
femtosecond pump-probe experiments where the damp- 
ing time is longer but also in low frequency Raman ex- 
periments where a better separation of scattering lines 
is possible. Second, as was shown by the pump-pump- 
probe experiments in Del Fatti et aLjl6| it is possible to 



control the acoustic motion of silver nanoparticles. With 
higher Q factor, the first pump pulse would "launch" the 
motion of every nanoparticles but the next pump pulse 
could selectively enhanced the vibration amplitude of a 
smaller subset of the nanoparticles. 

Solid matrix-embedded nanoparticles always have Q 
factor lower than free ones. However, free nanoparticles 
suffer from many drawbacks. For example, surface ad- 
sorption can affect their properties. Also it is often tech- 
nologically easier to deal with matrix embedded nanopar- 
ticles. Where solid embedding is technologically desir- 
able, our results provide guidance about how the matrix 
can be optimally chosen. 

Finally, the sphere radius is irrelevant in this paper 
{Q does not depend on i?). So there are some potential 
applications for other kinds of systems. Such damped 
vibrational modes of a virus in water may be relevant to 
its mechanical response to ultrasound. 17] 

All of the foregoing discussion applies only to 
spheroidal vibrations with ^ = 0. A similar analysis can 
be carried out for torsional modes, but a separate expres- 
sion for acoustic impedance is required for each value 
of £. It is not apparent how to make the extension to 
spheroidal modes with £ ^ 0. Calculations made on CdS 
spheres embedded in silica in a previous work0 revealed 
lightly damped fundamental torsional modes which we 
could not understand at that time. The situation is even 
more dramatic for gold spheres inside diamond for which 
the lowest torsional modes have very high quality fac- 
tors (Q - 900 for (TOR,^ = l,n = 0), Q - 20000 for 
(TOR,^ = 2,71 = 0), . . . ) 
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